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This  paper  studies  the  dynamic  behavior  of  doubly  fed  induction  generators  (DFIGs)  modeled  with  the  stator 
flux  orientation  regarding  rotor-  and  grid-side  converter  control  strategies.  The  investigations  include  modal 
and  sensitivity  analysis,  and  using  them  to  identify  and  characterize  the  instability  mode,  stator  dynamics 
analysis,  and  time  domain  simulations.  The  paper  first  deals  with  the  modal  analysis  of  the  DFIG  for  different 
control  strategies  and  operating  conditions.  Then  the  dynamic  and  transient  performance  of  the  DFIG  under 
voltage  dips  and  large  disturbances  are  investigated,  and  the  effects  of  rotor  and  grid-filter  control  strategies  on 
the  DFIG  transient  behavior  are  examined.  Increasing  the  rotor  and  grid-filter  closed  loop  bandwidths,  rotor 
current  active  damping,  and  compensation  of  back-emf  voltages  can  improve  the  rotor  and  grid-filter 
dynamics,  decrease  tracking  error  and  limit  the  rotor  current  transients.  However  they  can  lead  to  poorly  stator 
modes  and  deteriorate  the  DFIG  transient  performance.  Thus,  appropriate  selection  of  rotor  controller 
parameters  has  an  important  effect  on  the  DFIG  transient  performance,  and  consequently  on  the  future 
interconnected  power  system. 
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1.  Introduction 

There  are  several  types  of  adjustable  speed  generators  used  in 
wind  turbines  [1],  The  DFIG  is  a  commonly  used  variable  speed 
generator  in  wind  turbine  applications.  In  Fig.  1,  the  DFIG-based  wind 
turbine  is  connected  to  the  infinite  bus  through  the  equivalent  grid 
impedance.  It  consists  of  a  wound  rotor  induction  generator  with 
back-to-back  voltage  source  converters  linking  the  rotor  to  the  grid. 
The  rotor-side  converter  is  connected  to  rotor  windings  and  is  used  to 
control  the  generator  speed  and  reactive  power,  while  the  grid-side 
converter  is  connected  to  the  grid  through  a  grid-side  filter  and  is  used 
to  control  dc  bus  voltage  and  reactive  power. 

Field  oriented  control  (FOC)  is  the  most  widely  used  controlling 
method  of  DFIG,  which  enables  decoupled  control  of  real  and 
reactive  power.  The  FOC  may  be  implemented  as  stator-flux 
orientation  control  [2],  or  air-gap  flux  orientation  [3], 

In  recent  years,  DFIG-based  wind  turbines  have  been  an 
integrated  part  of  distributed  generation  system.  Therefore,  they 
will  have  strong  effects  on  the  future  interconnected  power 
system.  In  other  words,  any  abnormality  associated  with  the  DFIG 
affects  the  power  grid  performance.  Taking  this  into  account,  with 
increased  penetration  of  wind  power,  the  behavior  of  DFIG-based 
wind  turbines  under  disturbances  becomes  more  important. 

This  work  deals  with  the  dynamic  behavior  investigation  of 
doubly  fed  induction  generators  with  regard  to  rotor  and  grid- 
filter  control  strategies. 

The  dynamic  and  transient  studies  of  the  DFIG-based  wind 
turbines  have  been  carried  out  by  various  authors  that  can  be  cate¬ 
gorized  in  several  groups,  such  as  impact  of  DFIG  on  power  system 
dynamics,  grid  fault  response  of  DFIG,  low  voltage  ride-through  (LVRT) 
capability  enhancement  of  DFIG,  participation  of  DFIG  in  grid  fre¬ 
quency  regulation,  and  investigation  of  DFIG  dynamic  behavior. 

In  [4-10],  the  impact  of  DFIG  on  the  power  system  electro¬ 
mechanical  oscillations  has  been  studied.  In  [4-5],  small  signal 
stability  of  power  system  with  DFIG  is  investigated.  The  impact  of 
grid  connected  wind  farms  on  power  system  oscillations  is 
reviewed  in  [6],  In  [7,8],  the  DFIG  contributes  to  the  network 
oscillation  damping  with  using  evolutionary  techniques.  In  [9],  the 
DFIG  enhances  the  power  system  damping  via  network  auxiliary 
power  system  stabilizer  loop.  The  effect  of  DFIG  on  transient 
stability  margin  of  synchronous  generator  is  assessed  by  [10], 

Some  researches  have  been  carried  out  addressing  the  grid 
fault  response  and  short  circuit  current  of  DFIG  [11-16].  These 


papers  mainly  consider  the  short  circuit  behavior  of  the  DFIG,  and 
develop  theoretical  analyses  of  the  short  circuit  current. 

In  [17-20],  the  LVRT  capability  of  the  DFIG  has  been  investi¬ 
gated  in  order  to  limit  the  rotor  fault  current  and  dc-link  over 
voltage  within  their  acceptable  ranges  during  the  grid  voltage  dip. 
By  improving  the  LVRT  capability,  the  DFIG  remains  connected  to 
the  grid  during  and  after  clearing  the  fault,  and  contributes  to  the 
system  stability. 

Participation  of  the  DFIG  in  the  power  system  frequency  regula¬ 
tion  has  been  studied  in  [21-23].  Therefore,  DFIG  can  contribute  to 
the  grid  stability  enhancement  by  frequency  regulation. 

Various  researchers  have  dealt  with  the  dynamic  behavior  of 
the  DFIG  wind  turbines  [24-33],  In  [24],  dynamic  performance  of 
fixed  speed  wind  turbine  and  DFIG  is  compared  during  power 
system  disturbances.  Transient  stability  of  the  DFIG  due  to  short 
circuit  fault  is  studied  by  [25].  The  DFIG  torsional  oscillations  due 
to  wind  turbine  drive  train  dynamics  are  presented  in  [26-28],  In 
the  dynamic  studies  of  [29-30]  the  stator  dynamics  are  not  taken 
into  account,  and  in  [31-32]  the  dynamics  of  controllers  are  not 
considered.  Also,  in  [33]  there  is  no  discussion  about  the  critical 
modes  and  their  effects  on  system  dynamics. 

The  nature  of  transient  instability,  factors  influencing  it,  mode 
of  instability  and  the  effects  of  system  and  controller  parameters 
on  DFIG  dynamic  behavior  are  not  veiy  clear  in  time  domain  DFIG 
related  literatures.  Also,  little  work  has  been  carried  out  regarding 
poorly  damped  modes  of  DFIG,  and  their  effects  on  DFIG  dynamic 
performance. 

The  aim  of  this  paper  is  to  provide  insight  about  the  dynamic 
behavior  of  DFIG  based  wind  turbines.  Controller  parameters  of 
the  DFIG  have  an  important  impact  on  the  dynamic  performance 
of  the  wind  turbine.  Therefore  DFIG  controllers  will  play  a 
significant  role  in  the  future  power  system  with  increased  pene¬ 
tration  of  DFIG.  The  paper  considers  generic  PI  controllers  for 
regulation  of  rotor  and  grid-filter  currents,  and  investigates  the 
impacts  of  rotor  and  grid-filter  controller  parameters  on  the  DFIG 
dynamic  performance  under  voltage  dips. 

The  structure  of  the  paper  is  as  follows.  Following  this  intro¬ 
duction,  dynamic  model  of  the  DFIG  is  given.  Then,  modal  analysis 
of  the  system  is  presented.  Based  on  this  analysis,  poorly  damped 
modes  of  the  system  are  extracted  and  the  effects  of  these  modes 
on  the  DFIG  dynamic  performance  are  examined  both  analytically 
and  by  simulations.  Also,  transient  behavior  of  the  DFIG  with 
regard  to  different  control  parameters  is  examined. 


Converter  Converter 


Fig.  1.  DFIG  connected  to  the  infinite  bus. 
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List  of  symbols 

(Os  n  (Of 

(0 

synchronous  and  rotor  angular  frequency  (in  pu) 
angular  frequency  of  stator  flux  vector  (in  pu) 

Vrx,  voltage  amplitude  of  infinite  bus 

S  ,  (0  2 

rotor  slip  and  rotor  slip  frequency 

Vsdq  >  Wrdq  dq  components  of  stator  and  rotor  fluxs 

Ls,Lr 

stator  and  rotor  self  inductances 

irdq ,  isdq  igdq  dq  components  of  rotor,  stator  and  grid-filter  currents 

Ks  ,  Kr 

stator  and  rotor  resistances 

vrdq,  v^q  vadq  dq  components  of  rotor,  stator  and  grid-filter 

Lm 

magnetizing  inductance 

voltages 

s ,  r 

first  subscript  indicates  stator  and  rotor 

irdq _ rtf,  ivdq-ref  dci  components  of  rotor  and  grid-filter  refer- 

d  •  q 

second  subscript  indicates  direct  and  quadrature  axes 

ence  currents 

0 

subscript  corresponding  to  operating  point  quantities 

Pe  DFIG  electrical  real  power 

t ob  base  angular  frequency  (in  rad/s) 

g 

first  subscript  indicates  grid-side  filter 

2.  DFIG-based  wind  turbine  modeling  in  stator  flux  orientation 

The  purpose  of  this  section  is  to  present  the  dynamic  model  of 
single  machine  infinite  bus  (SMIB)  system  of  Fig.  1  in  synchronous 
reference  frame  with  the  stator  flux  orientation.  The  discussion 
considers  generic  PI  controllers  for  regulation  of  rotor  and  grid- 
filter  currents,  rotor  speed  and  reactive  power.  The  generalized 
machine  model  is  developed  based  on  the  following  conditions 
and  assumptions: 

(a)  Positive  direction  for  the  stator  and  rotor  currents  is  assumed 
into  the  generator  (see  Fig.  1). 

(b)  The  equations  are  derived  in  synchronous  reference  frame 
using  direct  (d)  and  quadrature  (q)  axes  representation. 

(c)  All  system  parameters  and  variables  are  in  per  unit  and 
referred  to  the  stator  side  of  DFIG. 


The  stator  and  rotor  voltages  and  fluxes,  electromechanical 
torque,  and  reactive  power  injected  to  the  grid  by  the  DFIG,  in  the 
reference  frame  rotating  at  angular  speed  of  co,  are  given  by  [17]: 


.  1  di//sda 

vsdq  —  Rs  hdq  +, JMW  sdq  +  ^  (f( 

(1) 

1  dwrda 

^rdq  Rrbdq  ~^J^2Wrdq  “1“  fa 

(2) 

1/7 sd  =  Lshd  +  Lmlrd  l// Sq  =  Lslsq  +  Lm  lrq 

(3) 

Wrd=  bd  “1"  Lr  ird  1/7  ,.g  =  Lm  iSq  +  Lr  irq 

(4) 

L e  =  -j —(Wsqbd  ~  W sdbq) 

(5) 

Qs  =  Vsd-hq  ~  Vsq-hd 

(6) 

2.1.  Rotor  modeling 

From  (2)  to  (4),  the  rotor  dynamics  is  described  in  terms  of 
rotor  current  and  stator  flux,  as  follows: 

J  /  ^  | 

~  =  ~Erirdq—j(tJ2Lrdrdq—Edq-\-Vrdq  (2) 

I2  n  \2 

where,  L'r  =  Lr  -  ,  R'r  =  Rr  +  f j-  J  Rs,  and 

Edq  =j)(Vsdq -jOJrWsdq-J^Wsdq^)  (8) 


In  order  to  decrease  tracking  error,  the  back-emf  voltages,  Ed 
and  Eq,  can  be  compensated  by  rotor  current  controllers  using  feed 
forward  terms.  Considering  the  rotor  controllers  to  be  PI, 
Kidq(s)  =  kp-idq  +  l<i-idq/s<  and  under  compensation  of  cross  cou¬ 
pling  terms,  the  d-q  rotor  voltage  could  be  stated  as 

^rdq(f)  =  Ep  idq-Ordq-ref(t)  ~  Irdq(O) 

+  idq  ■  jih'dq  ref(L)  —  irdq(t))dt  +  jO)2Lrirdq  +  kr  _  cam-Edq(t) 

(9) 


Considering  (9)  as  the  rotor  control  voltage,  the  rotor  d-q 
current  control  loops  can  be  described  by  Fig.  2(a).  In  this  figure, 
the  back-emf  voltages,  Ed  and  Eq,  are  represented  as  disturbances. 
In  (9)  and  Fig.  2(a),  kr_com  is  either  0  or  1.  fcr_com  =  1  means  that 
back-emf  voltages  are  compensated  by  rotor  current  controllers, 
and  kr-com  =  0  means  that  they  are  not  compensated.  Considering 
Fig.  2(a),  the  open  loop  bandwidth  of  current  control  in  per  unit  is 
as  =  R'r/L'r,  which  is  relatively  small.  By  pole  cancellation  of  the 
plant  with  the  zero  of  controller,  in  Fig.  2(a),  we  have 


kj  -  idq 
kp  -  idq 


(10) 


Selecting,  kp_idq  =  adqL'r/a>b,  and  by  taking  (10)  into  account, 
the  transfer  function  from  irdq^ref  to  iriq  will  be 


G 


cl  -  idq 


Idq(S) 
Idq  -  re/(S) 


&dq 

S-\-Otdq 


(ID 


adq  is  the  closed  loop  bandwidth  of  d—q  rotor  current  control, 
adq=(kp_idq/L'r)a>b.  Also,  the  transfer  function  from  Edq  to  irdq  is 
given  by: 

Irdq(S)  COb  S 

Edq(S)  L'r  (s+adq)(s+as) 


Normally,  as  is  small  and  therefore  the  rotor  current  response  to 
back-emf  voltages  is  relatively  slow.  To  overcome  this  problem,  the 
open  loop  bandwidth,  as,  can  actively  be  increased  by  means  of 
control,  this  is  called  active  damping.  The  current  control  loop  with 
active  damping  is  depicted  in  Fig.  2(b).  Considering  Fig.  2(b),  the 
improved  open  loop  bandwidth  will  be  as_active  =  (R'r+Rr_active/ 
L'r)a>b.  From  (7)  and  (9)  and  Fig.  2(b),  the  rotor  current  dynamics  can 
be  presented  as 


1,-  dird 
wb  dt 


LLrird~E(kr  —  com  1  )Ed  -F  kp  ,p(  ird  Ed)  T  ^5  Rr  activeEd 

(13) 


Ed  and  Eq  in  (8)  are  functions  of  stator  flux  and  stator  voltage.  These 
terms,  called  rotor  back-emf  voltages,  reflect  the  effects  of  stator 
dynamics  on  rotor  current  dynamics  and  have  an  important  role  on 
DFIG  transient  performance.  By  compensating  the  cross  coupling 
terms,  co2L'rirq  and  - co2L'rird,  by  d-q  rotor  current  controllers,  the  d 
and  q  rotor  current  control  loops  will  be  decoupled. 


Lr  dirq 
wb  dt 


-  -  Rrirq  +  ( kr 


1  )Eq  +  kp _  iq(irq _ ref  —  ifT/t  +  Xp  —  R, 


activebq 


(14) 


dx 5 
dt 


-  id  Ord  -  ref  bd  ) 


(15) 
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l _ — _ J 

d-axis  rotor  current  control  d-axis  rotor  current  dynamics 


I _ J 


q-axis  rotor  current  control  q-axis  rotor  current  dynamics 


Fig.  2.  (a)  Rotor  current  control  loop,  (b)  rotor  current  control  loop  with  active  damping. 


dxQ 

-jjr=  ki-iq(irq-ref  ~irq)  (16) 

Rr-  active  is  the  rotor  active  resistance  caused  by  active  damping  shown 
in  Fig.  2(b). 

2.2.  Grid-side  filter  modeling 

The  grid-side  filter,  as  shown  in  Fig.  1,  consists  of  inductance  Lg 
and  resistance  Rg,  and  its  dynamics  are  described  by: 

-  Jp  =  -  Rg  igdq  igdq  —  Vgdq  +  Vsdq  (17) 

where  vsdq ,  igdq  and  vgdq  are  d  —  q  components  of  the  generator 
terminal  voltage  and  grid-side  filter  current  and  voltage, 


respectively.  vg  is  provided  by  the  grid-side  converter.  By  com¬ 
pensating  the  cross  coupling  terms,  ajLgigq  and  -  coLgigd,  the  d  and 
q  grid-filter  current  control  loops  will  be  decoupled.  Fig.  3  depicts 
the  d—q  grid-filter  current  control  loops.  The  voltages  vsd  and  vsq, 
in  Fig.  3,  are  represented  as  disturbances  in  the  current  control 
loops.  In  order  to  decrease  tracking  error,  these  voltages  can  be 
compensated  using  feed  forward  terms. 

Considering  Fig.  3  and  assuming  the  grid-filter  current 
controller  to  be  PI,  I<Cdq(s)  =  kp_g+kj^g/s,  the  d—q  grid-filter  voltage, 
under  compensation  of  cross  coupling  terms,  could  be  stated  as 

Vgdq  (0  —  ~  kp-g(  igdq  -  re/(  t)  —  igdq  ( t)) 

i-g  J  (igdq -refit)  ~  igdq(i-))dt —jGdLgigdq+kg_comVsdq(t)  (18) 
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d-axis  grid  filter  current  control 


d-axis  grid  filter  current  dynamics 


K sq - 


K„.c 


IVv, 


kgJ?—C> 


■gq 


V 


-Vgq!  J 

1 

1 . 

1  *§<1 

- i - < 

—  s+R„ 

coh 

1 

l. 


Igd  ■ 


q-axis  grid  filter  current  control  q-axis  grid  filter  current  dynamics 

Fig.  3.  Grid-side  filter  current  control  loop. 


In  (18),  kg  _com  =  1  means  that  back-emf  voltages  are  compen¬ 
sated  by  the  rotor  current  controllers,  and  kg_com  =  0  means  that 
they  are  not  compensated.  From  (17)  and  (18),  the  grid-side  filter 
dynamics  can  be  presented  as 

=  —  Rgigd  “F  kp~g(igd~ref  ~  lgd)“F -*14  +  (1  —  kg_com)Vscj  (19) 

(Of,  dt  =  ~^g<l~k^P-g^gq-ref  —  igq)+x^5  +  ti—kg-c0m)Vsq  (20) 

— [jr=  ki-g(igd-ref  —  !gd)  (21) 

—^=l<i-g(igq-ref  —  ^gq)  (22) 

The  PI  controller  parameters  of  rotor  and  grid-filter  current 
controls,  in  (9)  and  (18),  are  selected  such  that  the  closed-loop 
bandwidths  of  the  control  loops,  in  Figs.  2  and  3,  are  identical. 


2.3.  Stator  modeling 


In  stator  flux  orientation,  y/s  =  y/slt  and  y/sq  =  0.  Then,  from  (1), 
(3)  and  (4),  the  stator  can  be  described  by  the  following  state  equations 
as  a  function  of  rotor  current,  stator  flux,  and  infinite  bus  voltage. 


1  Ls  +  he  tty/sd 
cob  Ls  dt 
LeU 
Ls 


Rs3~Re  Rs~kRer 

r  sd  - j - wn  ‘  id  ~  Iv-  igd 


(Oirq  -{-  LpCOigq  4- 


LqL, 


Ls 

,  J  dira 
a>b  dt 


LP  di. 


gd 


cob  dt 


-l-Vso  cos  y 


(23) 


dy 

-^  =  c0b(a>s-a>) 

and, 

,  _RS+Re,  .  p  hq+utldiS-^+v^sinr 
CO  —  z  Lm  lrq  —  Ke 


Vsd  ( 1  +  ij)  —  J^Lm  bd  +  Le  igd 


(24) 


(25) 


where,  cos  =  (L/a>b)(d0s/dt)  and  w  =  (1  / u>b)(dd / dt).  0S  and  0  are  the 
infinite  bus  voltage  angle  and  stator  flux  angle  in  stationary  reference 
frame,  respectively.  Also,  y  is  the  difference  between  0S  and  0. 


2.4.  Drive  train  model  and  speed  controller 


The  drive  train  comprises  turbine,  gear  box,  shafts  and  other 
mechanical  components  of  wind  turbine.  A  multi  mass  (usually 
two  mass)  model  must  be  used  for  dynamic/transient  studies  of 
wind  turbines  with  DFIC.  The  two  mass  model  is  given  by  [31], 

d(Dr  Te  +  ksp +D(a>t-  wr)  , 

-df= - 2Hr -  (26) 

da>[  _Tm-  ksp  -D(a>t-  wr) 

dt  ~  2 Hc  K  ’ 


dp 

dt 


COb(wt-COr) 


(28) 


where  wt  and  wr  are  the  turbine  and  generator  speeds  in  (pu), 
P  is  the  shaft  twist  angle  in  (rad),  Hr  and  Ht  are  the  inertia 
constants  of  turbine  and  generator  in  (s),  respectively,  ks  is  the 
shaft  stiffness  coefficient  in  (pu/elec.  rad),  D  is  the  damping 
coefficient  in  (pu),  Te  and  Tm  are  the  generator  electrical  torque  and 
the  turbine  mechanical  torque,  respectively,  in  (pu).  With  stator  flux 
orientation,  the  rotor  speed  is  controlled  by  the  q-components  of 
rotor  voltage  and  current,  vrq  and  irq  [1],  The  control  scheme  used 
for  speed  control  is  shown  in  Fig.  4.  In  this  figure,  aq  is  the  bandwidth 
of  the  q-axis  rotor  current  control  loop,  and  Tdis  =  ks0+Da)f  Employ¬ 
ing  a  PI  controller  for  the  speed  controller,  !<a,(s)  =  kp^  +  k/ai/s,  state 
equation  of  the  speed  controller  is 

i .  —  kfr,,(rur  _  o>i)  (29) 


The  q-axis  rotor  current  reference,  for  speed  control  and 
damping  of  rotor  oscillations  is 

Irq-ref  =  ~kpto{(l>r^ref  —  OJr)  —  xl  (30) 


2.5.  Reactive  power  control 

With  stator  flux  orientation,  terminal  voltage  and  reactive 
power  exchange  between  the  generator  and  the  grid  can  be 
controlled  by  the  d-components  of  rotor  voltage  and  current  [1], 
According  to  (6),  in  vector  control  with  stator  flux  orientation,  the 
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Fig.  4.  Speed  control  loop. 


Table  1 

Stator  power  factor  as  a  function  of  d-axis  rotor  current. 


ird0(pu)  0  0.0862  0.172  0.345  0.52  0.69  0.862 

Power  factor  0.94  lag  0.96  lag  0.984  lag  1  0.98  lead  0.92  lead  0.879  lead 


Vs 


Fig.  5.  Reactive  power  control  loop. 


reactive  power  injected  to  the  grid  by  the  stator,  Q.s,  can  be  written 
as 

QS  =j\,  Wsd(Lmhd  ~  VGd)  (3D 

As  can  be  seen,  there  is  direct  relation  between  d-axis  rotor 
current  and  stator  reactive  power  and  the  generator  power  factor. 
Table  1  shows  the  stator  power  factor  with  different  d-axis  rotor 
currents,  for  DFIG  system  of  Appendix  A. 

Thus,  the  d-axis  rotor  reference  current  is  determined  by 
reactive  power  controller,  as  shown  in  Fig.  5.  ad  in  Fig.  5  is  the 
bandwidth  of  the  d-axis  rotor  current  control  loop.  Using  the  PI 
controller  as  Kq(s)  =  kp-q  +  kj-q/s,  the  state  equation  of  the 
reactive  power  controller  is 

dx 

-^=fc;-a«L-re/-Qs)  (32) 

The  d-axis  rotor  current  reference,  for  reactive  power  control  is 

ird-ref=kp—Q(Q.s-ref  Qs  )“FX8  (33) 


3.  Small  signal  stability  and  modal  analysis 

Eqs.  (13)-(16),  (19)-(30)  and  (32), (33)  describe  the  dynamics 
of  turbine-generator  and  rotor  speed  and  reactive  power 
controllers.  The  dynamic  model  of  the  DFIG  wind  turbine  may 
be  given  by 

x‘  =f(x,z,  d ) 

0  =  g(x,z,d)  (34) 

where  x,  z  and  d  are  the  vectors  of  the  DFIG  state  variables, 
reference  inputs  and  exogenous  inputs,  respectively,  i.e., 

X  =  Wsd>  /.  hd,  irq,X5,X6,  igd,  igq,Xu,X\5,X7,Xg,  Wr,  9,  Wt]T, 

Z  —  [Yvr  ,£?/-  Q.S  -ref]  >  d=[VOQ,  Tm ] 


Table  2 

The  sample  DFIG  system  modes  and  dominant  state  variables. 


System  modes  State  variables  with  highest  participation  factors 


h,2 

-7.03  +J328.8 

Wsd 

r 

h 

-755.6 

>rd 

irq 

-735.6 

l'r, 

ird 

^5,6 

-4+J12.69 

COr 

P 

A7 

-4.54 

*8 

- 

^8,9 

-0.57+J1.56 

*7 

cot 

AlO 

-13.27 

*5 

- 

An 

-14.7 

*6 

- 

Al2,13 

-754 

igd 

lgq 

Al4,15 

-3.63 

XM 

*15 

Linearizing  and  rearranging  (34)  yields  the  linearized  model  of 
the  DFIG-based  wind  turbine  as  follows. 

Ax‘  =A.Ax+B.Au  (35) 

3.3.  Modal  analysis  without  compensation  of  rotor  and  grid-filter 
back-emf  voltages 

Modal  analysis  is  carried  out  on  the  SMIB  system  shown  in 
Fig.  1.  Appendix  A  gives  the  parameters  (in  per  unit)  of  the 
1.66  MVA,  575  V,  60  Hz  DFIG  used  for  the  study.  The  study  is  done 
under  operating  conditions  of  Appendix  B.  These  conditions 
correspond  to  unity  power  factor  at  the  stator  terminal.  The  PI 
controller  parameters  are  given  in  Appendix  C.  These  parameters 
correspond  to  the  rotor  and  grid-side  filter  current  control  band- 
widths  of  2  per  unit  (754  rad/s),  speed  control  loop  bandwidth  of 
4.4  rad/s  (0.7  Hz)  including  active  damping,  and  reactive  power 
control  loop  bandwidth  of  4.4  rad/s  (0.7  Hz).  Also,  the  grid 
impedance  is  assumed  to  be  Re+jXe  =  0.01  -fjO.l.  Table  2  shows 
the  system  modes  and  corresponding  state  variables  with  the 
highest  participation  factors  under  condition  that  rotor  back-emf 
voltages  are  not  compensated,  i.e.  with  kr-Com  =  0  in  (13)  and  (14). 
Using  the  participation  factors  [34],  the  contribution  degree  of 
each  state  variable  in  the  system  modes  and  the  physical  nature  of 
dynamic  modes  can  be  detected  (see  Appendix  D).  Considering 
Table  2,  the  following  key  points  are  found: 

1.  The  dynamics  of  DFIG  contain  poorly  damped  modes 
(-7.03  +J328.8)  with  a  corresponding  natural  frequency  close 
to  the  network  frequency.  Stator  variables,  y/si  and  y,  have  the 
highest  contribution  in  these  modes,  and  consequently  they 
represent  the  stator  modes.  As  will  be  shown  later,  these 
modes  have  significant  impact  on  transient  performance  of 
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DFIG,  and  under  special  operating  conditions  may  become 
unstable. 

2.  The  modes  A3  =  -755.6  and  14  =  -735.6  are  the  d—q  rotor 
current  modes,  and  ird  and  irq  have  the  highest  participation  in 
these  modes.  These  modes  are  very  fast  and  their  damping  is 
approximately  equal  to  the  rotor  current  control  bandwidth 
(2  pu  or  754  rad/s).  Thus,  the  larger  rotor  current  control  band¬ 
width,  the  larger  damping  of  rotor  current  modes. 

3.  The  modes  25,6  =  -4  +  jl2.69  are  the  electromechanical 
modes.  The  mechanical  variables  cor  and  /?  have  the  highest 
contribution  in  these  modes.  The  corresponding  natural  fre¬ 
quency  is  relatively  2  Hz. 

4.  The  real  mode  27  =  -4.54  is  associated  with  state  variable  x8, 
see  (32),  and  is  equal  to  bandwidth  of  reactive  power 
control  loop. 

5.  The  modes  28,9  =  —  0.56+jl.56  are  the  mechanical  modes 
associated  with  the  state  variables  x7  and  cot.  These  modes 
are  approximately  weak  and  are  dependent  on  the  speed 
control  bandwidth  and  damping. 

6.  The  modes  210  =  -13.27  and  2n  =  -14.7  are  the  rotor  elec¬ 
trical  modes  associated  with  state  variables  x5  and  x6  (see  (15) 
and  (16)).  These  modes  are  equal  to  the  d—q  rotor  current  open 
loop  bandwidth,  as  =  (R'r/L'r)wt,. 

7.  The  modes  Au.t3  =  -754  are  the  d-q  grid-side  filter  current 
modes  associated  with  state  variables  igd  and  igq.  These  modes 
are  very  fast  and  their  damping  is  approximately  equal  to  the 
grid-filter  current  control  bandwidth  (2  pu  or  754  rad/s). 

8.  The  modes  2i445  =  -3.63  are  the  grid-side  filter  modes  asso¬ 
ciated  with  state  variables  x14  and  x15  (see  (21)  and  (22)).  These 
modes  are  equal  to  the  grid-filter  current  open  loop  bandwidth, 
ag  =  (Rg/Lg)(i>h. 

As  stated  above,  the  stator  modes  are  weakly  damped,  e.g.  the 

damping  ratio  of  stator  modes  in  Table  2  is  £,  =  0.0213.  Usually  in 

Table  3 

The  sample  DFIG  system  modes  with  rotor  and  grid-filter  closed  loop  bandwidths 

of  4  pu. 


System  modes  State  variables  with  highest  participation  factor 


h,2 

-2.1  +J330.2 

Vsd 

r 

h 

- 1498.5 

bd 

irq 

X4 

-1513.6 

bq 

ird 

^5,6 

-3.9  +J12.9 

(Or 

P 

X7 

-4.4 

*8 

~ 

^8,9 

-0.6+J1.6 

*7 

cot 

^10 

-13.5 

*5 

- 

Xu 

-14.0 

*6 

- 

^12,13 

-950.2 

igd 

hq 

^14,15 

-3.7  +jo.3 

*14 

*15 

the  literature,  the  stator  dynamics  are  neglected.  However,  these 
modes  could  have  significant  effects  on  DFIG  transient  behavior. 

In  the  following  sensitivities  of  stator  modes  with  respect  to  rotor 
and  grid-filter  control  bandwidths,  rotor  active  damping,  rotor  and 
grid-filter  back-emf  voltages  and  DFIG  power  factor  are  examined. 


3.2.  Stator  mode  sensitivity  with  respect  to  rotor  and  grid-filter 
current  closed  loop  bandwidths 

Normally,  increasing  the  rotor  and  grid-filter  closed  loop 
bandwidths  can  improve  the  rotor  current  dynamics,  decrease 
tracking  error,  limit  the  rotor  current  transients  and  improve  the 
DFIG  low  voltage  ride-through  capability.  Participation  factors 
show  that  the  stator  modes  are  somewhat  related  to  the  rotor 
and  grid-filter  state  variables.  Table  3  shows  the  system  modes 
under  condition  that  rotor  and  grid-filter  closed  loop  bandwidths 
are  selected  as  4  pu.  Comparing  Tables  2  and  3  shows  that 
increasing  the  rotor  and  grid-filter  closed  loop  bandwidths  while 
keeping  all  other  parameters  constant,  decreases  the  damping  of 
stator  modes.  Fig.  6  shows  the  stator  modes  real  part  as  a  function 
of  rotor  and  grid-filter  closed  loop  bandwidths.  Thus,  varying  the 
rotor/grid-filter  control  bandwidths  causes  significant  changes  in 
damping  of  stator  modes. 

Also,  increasing  the  rotor  resistance  passively,  i.e.  increasing  R'r, 
or  actively,  i.e.  increasing  Rr- active  in  Fig.  2(b),  can  decrease  the 
stator  modes  damping.  Table  4  shows  the  system  modes  and 
dominant  state  variables  for  rotor  active  damping  term  of 
Rr  -  active  =  0.16.  Modal  analysis  shows  that  for  Rr- active  =  0.1 6  pu, 
the  stator  modes  become  -0.71+j310.  Thus,  including  active 
damping  in  the  rotor  current  control  loop  can  significantly 
decrease  the  stator  modes  damping. 


3.3.  Stator  mode  sensitivity  with  respect  to  d-axis  rotor  current 

As  it  is  clear  from  Table  1,  there  is  direct  relation  between 
d-axis  rotor  current  and  stator  reactive  power  and  consequently 
the  generator  power  factor.  Modal  analysis  shows  that  increasing 
the  d-axis  rotor  current  can  decrease  the  damping  of  stator  modes. 
Fig.  7  shows  the  real  part  of  stator  modes  as  a  function  of  d-axis 
rotor  current  for  two  different  values  of  rotor  and  grid-filter 
current  closed  loop  bandwidths.  Considering  Fig.  7,  it  is  clear  that 
increasing  the  d-axis  rotor  current,  while  having  large  bandwidths 
of  current  control  loops,  can  move  the  stator  modes  to  unstable 
state  and  deteriorates  the  DFIG  dynamic  performance.  In  other 
words,  moving  the  turbine-generator  to  leading  power  factor  can 
influence  the  DFIG  dynamic  behavior. 


Fig.  6.  Real  part  of  the  stator  modes  as  a  function  of  rotor  and  grid-filter  closed  loop  bandwidths. 
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3.4.  Modal  analysis  with  compensation  of  rotor  back-emf  voltages 

It  is  possible  to  include  a  feed  forward  compensating  term  in 
the  control  law  to  compensate  for  the  tracking  error  caused  by 
variations  in  the  rotor  back-emf  voltages.  Table  5  shows  the  system 


LpLrr 


.coirq  _  rej  -T  Lecoigl 


.  fef/7 

-ref  +  . 

Ls 


1  di 


rd  -  ref 


cob  dt 


Le  digd  _  rcj 
a>b  dt 


+  Voa  cos  y 


(36) 


dy_  (  ( Rs+Re/Ls)Lmirq_ref-Refq  ref  +  (LeLm/Ls)('l/a>b)(dirq_ref/dt)-(Le/CL>b)(digq„ref/dt)+VcosinY 

dt  y  l//s d(L  ~b(Le /Ls))  —  (Le/Ls)Lmird ^ref -b Leigd _ref 


modes  and  dominant  state  variables  after  compensating  of  rotor 
and  grid-filter  back-emf  voltages.  Comparing  Tables  2  and  5  shows 
that  compensation  of  back-emf  voltages  can  significantly  decrease 
the  damping  of  stator  modes,  but  the  other  system  modes  remain 
constant.  Therefore,  although  compensation  of  back-emf  voltages 
can  improve  the  DFIG  low  voltage  ride-through  capability  and 
limit  the  peak  value  of  rotor  current  transients,  but  it  may  decrease 
the  stator  modes  damping. 


The  equilibrium  points  of  these  equations  are  y0  and  y/sd0  rs  1. 
The  linearized  dynamic  model  of  (36)  and  (37)  around  the 
operating  point  is  obtained  as 


'^v'si 

"an 

di2 

"A^sd" 

Ajr* 

=  CUb 

a2i 

a22 

A  y 

in  which  au  to  a22  are  given  in  Appendix  E.  According  to  Appendix 
E,  the  exact  characteristics  polynomial  is  given  by 


4.  Analysis  of  stator  flux  dynamics 


If  the  current  control  loop  of  the  rotor  and  grid-filter  is  fast 
enough,  it  is  possible  to  decouple  the  flux  dynamics  from  other 
system  dynamics  and  set  the  rotor  and  grid-filter  currents  to  their 
reference  values.  This  means  that  to  study  the  stator  flux  dynamics, 
it  is  sufficient  to  consider  (36)  and  (37). 


1  Ls  +Le  dlj/sd 
cob  Ls  dt 


ks+kp  ks+ke  .  . 

I  v>sd  '  j  Gn -bd-ref  —  'v? >gd  -  ref 


2  4-  a>b 


/ks  +  ke  ^(k5  +  ke/Ls)i//sdo  —  (ks  +  ke/Ls)Lmirdo  +  ke/g^o\ 
\  Ls  LP  */Gdo(l  ~b(Le/Ls))~  (LeLm/Ls)irdo-\-Leigdo  / 


7ks+kg  VcodQ 

V  Ls  +  Le  D 


(39) 


where  v^q  =  cos  y0=  cos  y0,  v^q0  =  Vx  sin  y0,  and  D  =  y/Sd0 
(L+(Le/Ls))-(LeLm/Ls)ird0+Leigd0.  For  typical  operating  conditions, 
y0  is  close  to  {n/2),  Lm  ^  Ls  and  Le  <4LS.  Therefore,  the  simplified 


Table  4 

The  sample  DFIG  system  modes  with  rotor  active  damping  term,  Rr_aaivc  =  0.16. 


System  modes  State  variables  with  highest  participation  factors 


h,2 

-0.71  +J309.8 

Wsd 

Y 

h 

-755.98 

bd 

bq 

a4 

-730.8 

irq 

bd 

>*5,6 

-3.9  ±j\  3.05 

(Or 

P 

Ai 

-4.43 

*8 

- 

^8,9 

—  0.58  +jl.58 

*7 

cot 

>*10,11 

-216.13 +J21.55 

*5 

*6 

>*12,13 

—  474  +J131.3 

hd 

lgq 

>*14,15 

—  3.64  ±j0.53 

X14 

*15 

Table  5 

The  sample  DFIG  system  modes  with  compensation  of  rotor  and  grid-filter  back- 
emf  voltages. 


System  modes  State  variables  with  highest  participation  factors 


>*1,2 

-0.41  ±j376 

Vsd 

Y 

A3 

-  758.4 

bd 

- 

-  745.1 

bq 

- 

A5,6 

-3.89  +J13.11 

(Or 

P 

Ai 

-4.37 

*8 

- 

>*8,9 

—  0.57  +jl.56 

*7 

(Of 

>*10,11 

- 13.41 

*5 

*6 

>*12,13 

-753.98 

igd 

lgq 

>*14,15 

-3.77 

X14 

*15 

Fig.  7.  Stator  modes  real  part  with  respect  to  ira  under  different  rotor  and  grid  filter  control  bandwidths. 


rotor  voltage  magnitude  (pu)  Q_  rotor  current  magnitude  (pu)  O  Elec.  Torque  (pu)  O"  Stator  flux  magnitude  (pu) 
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Time  (sec) 


Time  (sec) 


Fig.  8.  Transient  behavior  of  DFIG  with  two  different  values  of  rotor  and  grid  filter  control  band  widths,  in  pu,  (a)  stator  flux  magnitude,  (b)  electrical  torque,  (c)  rotor  current 
magnitude  (d)  rotor  voltage  magnitude. 
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characteristic  polynomial  is  approximately  found  as 


/P  3-cOb 


(Rs+Re 

\Ls-\-Le 


(Rs Re  /  Ls)y/ sd  o  —  (Rs  +  Re/Ls)Lmird0  +  Reigdo\ 

WsdO^  ^  +  (fe  /  fs))  —  (  Wm  /  f s)*rdO  +  Le  lgdO  / 


+«g  =  o 
(40) 


Considering  (40),  it  is  clear  that  the  stator  mode  natural 
frequency  is  close  to  the  line  frequency,  i.e.,  1  pu.  For  maintaining 
stability  of  the  (38),  it  is  required  that 


'rd0  (1  +Ls  +  Le)  Lmlgd0  (is+Le  +  Rs+^) 


fm 


(41) 


If  inequality  (41 )  is  satisfied,  the  nonlinear  stator  dynamics  is 
locally  asymptotically  stable.  In  normal  conditions,  the  grid-side 
converter  operates  at  unity  power  factor,  and  consequently  d-axis 
grid  filter  current,  igd0,  exchanging  reactive  power  with  the  grid  is 
equal  to  zero.  Hence,  by  setting  igdo=0  in  (41),  we  have 


bdo 


2t/Jsd0 

lm 


(42) 


Considering  (42),  it  is  clear  that  stator  dynamics  is  dependent 
on  operating  condition,  ird0,  and  correspondingly,  stator  power 
factor.  Also,  it  depends  on  the  network  and  generator  parameters 
such  as  Ie  and  Ls. 


5.  Time  domain  simulations 

In  the  previous  section,  modal  analysis  of  the  DFIG  system  was 
presented  and  the  effects  of  different  rotor  and  grid-filter  control 
strategies  on  the  poorly  stator  modes  were  investigated.  The 
approximation  of  the  DFIG  dynamic  and  transient  characteristics 
can  be  obtained  from  the  eigenvalues  of  the  linearized  DFIG  model 
presented  in  Section  3.  In  this  section  dynamic  and  transient 
behavior  of  the  DFIG  system  is  studied.  The  simulations  are  carried 
out  on  the  SMIB  system  shown  in  Fig.  1.  The  system  and  controller 
parameters  of  the  DFIG  used  in  these  studies  are  shown  in 
Appendix  A  and  C.  The  DFIG  operating  slip  and  electrical  power 
before  the  fault  are  s0  =  -0.1  and  Pe0  =  0.5  pu,  respectively,  and 
the  generator  is  operated  at  unity  power  factor.  The  grid  impe¬ 
dance  is  assumed  to  be  Re+jXe  =  0.05+j0.5. 

5.1.  Effects  of  rotor  and  grid-filter  closed  loop  bandwidths 

To  find  out  the  effects  of  the  rotor  and  grid-filter  closed  loop 
bandwidths  on  the  DFIG  transient  behavior,  a  40%  voltage  dip  with 
duration  of  300  ms  is  imposed  on  the  DFIG  terminal  at  t=  0.9  s. 
Fig.  8  shows  the  transient  response  of  the  DFIG  with  two  different 
values  of  rotor  and  grid-filter  closed-loop  bandwidths.  Considering 
Fig.  8,  increasing  the  closed  loop  bandwidth  from  1  pu  to  4  pu 
limits  the  peak  value  of  rotor  current,  electrical  torque  and  stator 
flux,  but  it  results  in  oscillatory  transient  response  with  larger 
settling  time.  This  is  because  as  stated  in  Section  3,  the  stator  flux 
and  thus  the  DFIG  transient  response  is  less  damped  with  higher 
rotor  and  grid-filter  closed  loop  bandwidths.  Moreover,  it  is  clear 
from  Fig.  8  that  the  larger  closed  loop  bandwidth,  the  larger 
control  efforts  and  rotor  voltages  (vrd  and  vrq).  Therefore,  a  trade 
off  is  needed  between  the  transient  response  settling  time,  peak 
value  of  transient  responses  and  control  effort,  in  selection  of 
appropriate  rotor  and  grid-filter  closed  loop  bandwidths. 

5.2.  Effects  of  rotor  active  damping 

As  already  mentioned  in  Section  3.2,  including  active  damping 
in  the  rotor  current  control  loop,  the  term  Rr_ active  in  Fig.  2(b)  can 
significantly  weaken  the  stator  dynamics.  To  demonstrate  the 
impact  of  rotor  active  damping  on  the  DFIG  transient  behavior,  a 


50%  voltage  dip  with  duration  of  300  ms  is  imposed  on  the  DFIG 
terminal  at  t=0.9  s.  Fig.  9  shows  the  transient  response  of  the  DFIG 
with  and  without  rotor  current  active  damping  for  Rr-active  =  0.2 
and  Rr-active  =  0.  It  is  clear  from  Fig.  9  that  the  DFIG  stability 
margin  in  the  case  with  the  rotor  active  damping  is  lower  than 
that  of  the  case  without  active  damping. 

Therefore,  the  settling  time  and  oscillations  of  transient 
responses  in  the  case  with  Rr-active  =  0.2  are  much  higher  than 
that  of  the  case  with  Rr  -  active  =  0.  This  is  because  the  stator  flux 
and  thus  the  DFIG  transient  response  is  less  damped  with  higher 
rotor  resistance.  Hence,  increasing  the  rotor  active  damping 
cannot  improve  the  DFIG  low  voltage  ride-through  capability.  This 
in  turn  can  deteriorate  the  DFIG  transient  performance,  i.e.  larger 
oscillations  in  stator  flux,  electrical  torque,  terminal  voltage,  and 
rotor  currents. 

5.3.  Effects  of  rotor  back-emf  voltages 

To  find  out  the  effects  of  compensation  of  rotor  back-emf 
voltages  on  the  DFIG  transient  behavior,  a  simulation  is  carried 
out  for  a  50%  voltage  dip  with  duration  of  300  ms.  In  this 
simulation  (Fig.  10)  two  states  are  considered.  In  the  first  state 
the  back-emf  voltages  are  not  compensated,  and  in  the  second 
state  they  are  compensated.  It  is  clear  from  Fig.  10  that  compensa¬ 
tion  of  these  voltages  limits  the  peak  value  of  the  rotor  current 
transients,  but  it  decreases  the  DFIG  stability  margin.  Therefore,  by 
compensating  of  back-emf  voltages,  the  settling  time  and  oscilla¬ 
tions  of  stator  flux  and  electrical  torque  are  increased.  Moreover, 
this  state  requires  larger  rotor  voltages  and  control  efforts. 

5.4.  Transient  stability  study 

As  stated  in  Section  4  and  considering  (42),  the  stability  of 
stator  dynamics  is  closely  related  to  the  network  reactance,  Le.  To 
demonstrate  the  prominent  role  of  stator  dynamics  on  the  DFIG 
transient  performance,  a  70%  voltage  dip  with  duration  of  300  ms 
is  imposed  on  the  infinite  bus  at  t=0.9  s.  In  this  simulation,  Fig.  11, 
two  different  network  impedances  are  considered.  In  the  first  case 
Re  +jXe  =  0.05+j0.3pu,  and  in  the  second  case,  Re+jXe  = 
0.05+j0.5pu.  In  both  cases,  the  slip  and  real  power  of  the  system 
at  the  equilibrium  point  are  -0.1  and  0.5  pu. 

It  is  clear  from  Fig.  11  that  in  case  2  with  higher  network 
reactance,  the  stator  flux  at  the  time  of  clearing  the  fault  is  out  of 
the  stator  domain  of  attraction,  and  consequently,  stator  state 
trajectories  do  not  reach  their  stable  post  fault  operating  point  and 
transient  behavior  of  the  DFIG  is  unstable.  The  final  outcome  of 
transient  instability  is  large  oscillations  in  the  terminal  voltage, 
electrical  torque  and  rotor  speed  of  DFIG. 

In  wind  turbine-generators,  the  nature  of  instability  is  different 
from  the  rotor  angle  instability  of  conventional  synchronous 
generators.  In  the  wind  turbines  with  the  DFIG  the  generator 
speed  range  is  approximately  +  30%  around  the  synchronous 
speed.  The  upper  limit  of  the  generator  speed  is  determined  by  the 
back-to-back  converter  capacity  rating.  If  the  generator  speed  after 
clearing  the  fault  is  higher  than  the  limit,  the  converter  cannot 
handle  the  slip  power  and  the  generator  may  become  unstable. 
For  the  DFIG  system  studied  in  this  paper  the  operating  speed 
before  the  fault  is  1.1  pu  and  the  generator  speed  at  the  moment  of 
clearing  the  fault  (at  t=  1.2  s)  is  1.11  pu.  Thus,  the  generator  slip 
after  the  fault  is  within  the  allowable  range  and  the  back-to-back 
converter  is  able  to  handle  the  slip  power.  Moreover,  growth  of  the 
generator  speed  during  the  fault  is  relatively  low.  Therefore,  the 
system  does  not  face  angle/speed  instability.  In  the  second  case,  in 
spite  of  admissible  slip,  the  system  experiences  instability  that 
roots  in  the  stator  dynamics  and  exhibits  as  transient  voltage 
instability. 
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Fig.  9.  Transient  behavior  of  DFIG  with  and  without  rotor  active  damping,  Rr-active  =  0  and  Rr_ active  =  0.2,  in  pu,  (a)  stator  flux  magnitude,  (b)  electrical  torque,  (c)  rotor 
current  magnitude,  (d)  rotor  voltage  magnitude. 
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Fig.  10.  Transient  behavior  of  DFIG  with  and  without  compensation  of  back-emf  voltages,  in  pu,  (a)  stator  flux,  (b)  electrical  torque,  (c)  rotor  current  magnitude,  (d)  rotor 
voltage  magnitude. 


864 


M.  Rahimi  /  Renewable  and  Sustainable  Energy  Reviews  37  (2014)  852-866 


0.8  1  1,2  1.4  1.6  1.8  2  2.2 

Time  (sec) 


Fig.  11.  DFIG  transient  response  under  300  ms  voltage  dip,  with  two  different  values  of  network  reactance  {Re+jXe  =  0.05  +j0.3pu  and  Re  +jXe  =  0.05+j0.5pu),  in  pu, 
(a)  stator  flux  magnitude  ( if/sd ),  (b)  terminal  voltage  magnitude,  (c)  electrical  torque. 


6.  Conclusion 

This  paper  studied  the  dynamic  and  transient  behavior  of  wind 
turbine  with  DFIG  with  regard  to  rotor  and  grid-filter  control 
strategies.  The  investigations  include  modal  and  sensitivity  analysis, 
stator  dynamics  analysis,  and  time  domain  simulations.  The  discussion 
considered  generic  PI  controllers  for  the  regulation  of  rotor  and  grid- 
filter  currents,  rotor  speed  and  reactive  power.  It  was  shown  that  the 
dynamics  of  the  DFIG  contain  poorly  damped  modes  corresponding  to 
stator  modes.  Also,  it  was  found  that  the  rotor  back-emf  voltages, 
closed  loop  bandwidths  of  rotor  and  grid-filter  currents,  and  rotor 
current  active  damping  have  high  impacts  on  the  stator  modes  and 
thus  on  the  DFIG  dynamic  performance.  The  approximation  of  the 
DFIG  dynamic  and  transient  behavior  can  be  obtained  from  the  modal 
analysis  and  eigenvalues  of  the  linearized  DFIG  model. 


The  paper  results  show  that  the  dynamic  behavior  of  the  DFIG 
under  voltage  dips  is  strongly  affected  by  the  stator  dynamics. 
Increasing  the  rotor  and  grid-filter  closed  loop  bandwidths  and 
compensation  of  back-emf  voltages  can  improve  the  rotor  and 
grid-filter  dynamics,  decrease  tracking  error  and  limit  the  rotor 
current  transients.  At  the  same  time,  they  can  lead  to  poorly  stator 
modes  and  deteriorate  the  DFIG  transient  performance.  Also, 
including  active  damping  term  in  the  rotor  current  control  loop 
can  decrease  the  DFIG  stability  margin  with  less  damped  transient 
response.  Therefore,  the  dynamic  performance  of  the  DFIG  is  to  a 
large  extent  dependent  on  rotor  controller  parameters. 

In  the  near  future  wind  farms  will  be  vital  and  important 
sections  of  the  power  system  generation.  Flence,  control  strategies 
of  the  DFIG  and  consequently  wind  farm  will  influence  on  the 
future  power  system  operation.  Also,  they  will  play  a  significant 
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role  in  order  to  guarantee  stable  and  secure  operation  in  the  near 
future  with  increased  wind  power  penetration. 


Appendix  A.  Parameters  of  the  1.66  MVA,  575  V,  60  Hz, 
DFIG-based  wind  turbine: 

Vbase  =  575  V,  Sbase  =  1 .66  MVA 

fbase  =  60  Hz,  cob  =  2nf  base  =  377  rad/s,  Rs  =  0.00706  pu 
Rr  =  0.005  pu,  Ls  =  3.07  pu,  Lr  =  3.056  pu,  Lm  =  2.9  pu 
Lg  =  0.3  pu,  Rg  =  0.003  pu 
Hr  =  0.75  s,  Ht  =  4.3  s,  D=1.2pu 
ks  =  0.6  pu/elec.rad,  oos  =  1  pu 


Appendix  B.  Operating  conditions  used  for  modal  analysis 
in  Section  3.1: 

Peo  =  0.9  pu,  o  =  —0.21,  vsd0  =  0, 
vS9o  =  1  pu  i//sd0  =  1 .005  pu,  i//sq0  =  0, 
ir do  =  0.345  pu,  irqo  =  0.7874, 


Appendix  C.  Controller  parameters  used  for  modal  analysis  in 

Section  3.1: 

kj]  j[/q  —  0.633  jciq  =  8.5,  kp  g  =  0.6 
kj_g  =  2.26  kpa,  =  6.98  kim  =  0.04656 
fcp_Q  =  0.01235  k,_Q  =  4.656 


Appendix  D.  The  participation  factors  of  mode  i  are  given  by 

Pi  =  [PliP2i--Pni]r;  Pki  =  ^  =  Wik4>ki 

where  akj  is  the  element  of  jth  column  and  kth  row  of  state  matrix 
A.  pki  is  a  measure  of  the  relative  participation  of  the  kth  state 
variable  in  the  ith  mode,  and  vice  versa  [34], 

i//ik  is  the  kth  element  of  the  ith  left-eigenvector  of  state  matrix 
A,  and  rpki  is  the  kth  element  of  the  ith  right-eigenvector  of  state 
matrix  A. 


Appendix  E.  The  elements  cin.  au,  a2 1  and  a22  of  the  state 
matrix  in  Eq.  (38): 


an 


Rs  +  Pp  Le 

Ls~kLe  D 


a  12  = 
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